Inverse Spin-Galvanic Effect in a Topological-Insulator/Ferromagnet Interface 
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When a ferromagnet is deposited on the surface of a topological insulator (TI) the topologically 
protected surface state develops a gap and becomes a 2-dimensional quantum Hall liquid. We 
demonstrate that the Hall current in such a liquid, induced by an external electric field, can have a 
large effect on the magnetization dynamics of the ferromagnet by changing the effective anisotropy 
field. This change is dissipationless and may be substantial even in weakly spin-orbit coupled 
ferromagnets. We study the possibility of dissipationless current-induced magnetization reversal in 
monolayer-thin, insulating ferromagnets with a soft perpendicular anisotropy and discuss possible 
applications of this effect. 

PACS numbers: 



Introduction. — Understanding the electric-field con- 
trol of magnetization and harnessing its technological 
potential are amongst the most important objectives of 
spintronics. Current-induced spin torques can reverse 
the magnetization of conducting ferromagnets and move 
magnetic domain walls [![. However, the Joule heating 
generated by transport currents remains a handicap from 
a practical viewpoint. An electric field can also reorient 
the magnetization of insulating compounds with broken 
inversion symmetry via the magnetoelectric coupling Q . 
While they overcome the issue with Joule heating, these 
multifcrroic materials are fewer and more difficult to en- 
gineer than common metallic ferromagnets. Recently, 
a novel magnetoelectric effect has been discovered |3[ 
in topological insulators that are coated with ferromag- 
netic films. TIs are bulk insulators with an anomalous 
band structure that supports topologically robust gapless 
states at the surfaces [4|. These materials are predicted 
to display a variety of unconventional spintronics effects 
Q . One unique feature is the universal quantized topo- 
logical magnetoelectric effect described by 




(b) 



L FM iA rj 

1^4 



FIG. 1: Corbino-disk-shaped TI coated with an ultrathin fer- 
romagnet. (a) Top view: in absence of electric fields, the mag- 
netization of the ferromagnet points outside the page (dotted 
circles). When a voltage difference is applied between the in- 
ner and outer circles, a dissipationless Hall current flows at 
the interface between the two materials (solid arrows). This 
current magnetizes the surface states of the TI (inverse spin- 
galvanic effect) along the radial direction (dashed arrows), 
thus exerting a spin torque on the magnetization of the fer- 
romagnet. (b) Cross sectional view: the shaded region is the 
TI, whereas the unshaded region is the ferromagnetic film. 
Hfm is the anisotropy field in electric equilibrium. Hcs.i is 
a topological magnetic field proportional to (and parallel to) 
the applied electric field. 
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Here M top is the induced magnetization, C\ is a half- 
integer topological invariant that depends solely on the 
sign of the time-reversal-symmetry-breaking perturba- 
tion, E is the applied electric field and C\e 2 /2tt = o~h is 
the Hall conductance (H = 1 throughout). Unfortunately, 
the prospects for manipulating the magnetization of real 
ferromagnets via Eq. (JTJ arc limited because below the 
threshold Hall current density (Jh < lA/m, see Ref. @) 
the topological magnetic field B top = (J,oM top < 10 -6 T is 
very small compared to typical coercive fields (> 0.01T) 
in a ferromagnet. 

In this Letter we identify a new contribution to the 
topological magnetoelectric effect, which stems from the 
current-induced spin polarization of the TI surface states. 
Unlike Eq. (JTJ, this effect depends on material parame- 
ters and is not related to Ampere's law; instead it is the 



topological counterpart of the inverse spin-galvanic effect 
found in conducting materials 0] ■ In ultrathin (thickness 
< lnm) ferromagnetic insulators deposited on a surface 
of TI (Fig. [!} the topological inverse spin-galvanic effect 
leads to qualitatively stronger spin-torques than Eq. 
thus opening the avenue for current-induced control of 
magnetization without Joule heating. 

Functional integral formalism. — We begin by review- 
ing the equation of motion for the magnetization M = 
MCl of a classical ferromagnet (in units of 1 /volume). 
At low energies the magnitude M is approximately con- 
stant and the only dynamical variable is the direction 
O = (Cl x ,Q y ,Q, z ). The time dependence of may be 
determined using the functional integral approach Q|, 
which is built on the partition function 



Z = Z Dfl(x,t)e 
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Zq is the partition function corresponding to the equilib- 
rium magnetic configuration SI = Sl cq . Sfm = Sb — £ 
is the action for small (quadratic) spin fluctuations, 

where Sb = M J dxdtfl ■ (fi eq x f2) is the Berry phase 
and £[ti] = J dxdtClx -1 ^ is the micromagnetic energy 
functional. \ is the spin-spin response function. The 
semiclassical equation of motion can be derived from 
6S FM /SCl = 0, 
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A gradient expansion Q of \ yields the venerable 
Landau-Lifshitz-Gilbcrt-Slonczewski equation for mag- 
netization dynamics in the presence of damping and 
transport currents, 



After rotating the spins by an angle ir/2 around z, Eq. 
J7|) may be rewritten as Sti = / d 2 xdtip[do — fi — 1f\ip 
with 

H = vft x (tt x — ea x ) + VFT y (ir y — ea y ) — Q z At z , (8) 

where ip is the rotated fermion field. In this new basis, 
a = A/(eup)(f2 x z) appears as an additional contri- 
bution to the effective vector potential. fi z A acts as 
a mass term. These massive Dirac fermions may be 
integrated out in the standard manner (Toj . whereby 
Z = J £>Sl(x,£)e~ Srff ^]. To second order in SI the ef- 
fective action is S c s ~ Sfm + Scs + Seb, where 



Scs = ^C x J d 2 xdte^ x A^A> 



(9) 



SI = Sl oq xH — aSl cq xO — v s -VSl — /3Sl oq xv s -VS7- 



(4) 



H is an effective magnetic field (in energy units) that 
includes the anisotropy field, the exchange field as well 
as external magnetic fields. H determines the easy axis 
along which the magnetization of a single-domain fer- 
romagnct points in equilibrium. v s is the adiabatic spin 
transfer velocity and is proportional to the transport cur- 
rent, a and (3 characterize dissipativc processes in which 
energy is transferred from magnetic to non-magnetic (e.g. 
lattice) degrees of freedom. 

Topological effective magnetic field. — We now address 
the magnetization dynamics of an insulating ferromagnet 
sitting on top of a TI. The low-energy effective Hamilto- 
nian for the surface states of the TI is 0, 0] 



H = vpT ■ (it x z) — At • SI, 



(5) 



where vp is the Fermi velocity, t 1 (i <= {x, y, z}) are 
Pauli matrices denoting real spin of the surface states, 
7r = — zV — eA, A is the electromagnetic vector poten- 
tial, z is the unit vector normal to the interface between 
the TI and the ferromagnet and A is the exchange cou- 
pling between the surface states and the local moments 
of the ferromagnet (A > for ferromagnetic coupling). 
We consider a ferromagnet with perpendicular anisotropy 
(Sl oq = z) so that in equilibrium a gap opens in the en- 
ergy spectrum of the surface states. 

The partition function for this composite system is 

Z = Zv J Df2(x, i) e - s ™["] j D 2 *(x, t)e~ s ^^' h \ 

(6) 

where 5*fm is the ferromagnetic action discussed above 
and 
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d 2 xdm [do-fi-H]® 



(7) 



is the action for the surface states. ^ is a fermionic 
spinor, do = dt — eA$, fi is the chemical potential (lo- 
cated in the gap) and Aq is the electrostatic potential. 



A = (Aq, A x + a x , A y + a y ) is the effective vector poten- 
tial and \i = t,x, y. The Chern-Simons action ((9|) arises 
in (2+1) dimensional systems with broken time reversal 
symmetry and nontrivial topology. The topology of the 
band structure is encoded in the TKNN ll| invariant C\ . 
For fermions described by a single Dirac Hamiltonian (J8]) 



we have 12 



Ci = --sgn(S\A). 



(10) 



iSeb is quadratic in spatial and temporal derivatives of 
Ap, and encodes the ordinary dielectric/diamagnetic re- 
sponse of the insulator. Herein we focus on ScSi which 
is first order in the derivatives of _4 M and thus domi- 
nates over 5eb at long length and time scales. It also 
produces the effective magnetic field that underlies the 
inverse spin-galvanic effect which is central to this study. 

The semiclassical magnetization dynamics follows from 

ss cti /sn = o, 



SI = Sl cq x (H 
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where Hfm = — 6Sfm/(M8Ci) is the effective magnetic 
field that collects the anisotropy /exchange fields of the 
isolated ferromagnet and 



H 



1 <55, 



_^lA[ e + A (f xsV 

M 2 d SCI M 2 u evF L evF 

(12) 

is an additional (topological) contribution to the mag- 
netic field that results from the exchange coupling be- 
tween the ferromagnet and the TI. M 2 d is the areal mag- 
netization at the interface (in units of 1/area). Hcs de- 
pends on material parameters (vp, A, M 2 d) and is pro- 
portional to the Hall conductivity an = C\e 2 /2ir. Be- 
cause the exchange coupling between the surface states 
and the localized moments of the ferromagnet is local in 
space, the influence of Hcs weakens as the thickness of 
the ferromagnetic film increases. 
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FIG. 2: Feynman diagrams for (a) the electric-field-induced 
magnetization (inverse spin-galvanic effect), (b) the xy com- 
ponent of the spin-spin response function, (c) the xy com- 
ponent of the spin-spin response function in presence of an 
electric current (it yields the adiabatic spin transfer torque 
v s ■ VO). The solid straight lines are propagators for mas- 
sive Dirac quasiparticles (quasiholes) . The solid wavy lines 
are magnons that couple to the spin operator and the dashed 
straight lines are photons that couple to the velocity operator. 



Hcs,i = -A/(et)j7M 2 D)Tf/E can be interpreted as an 
electric-field induced change of magnetic anisotropy. The 
underlying cause of this effect is that the electric field 
spin-polarizes the surface states along a direction (E/E) 
which is misaligned with the equilibrium easy axis (z). 
We illustrate this point by computing the magnetization 
induced by a static and uniform electric field: 



(13) 



where 



Xm.e 



(14) 

is the linear magnetoelectric response function (Fig. [2^). 
n, n' are the band indices of the surface states, i?k,n 
are the band energies, /k.„ are the Fermi distributions, 
T n n' = ( n > k|r l |n', k), and A is the area of the interface. 
From Eq. (J5J), the velocity operator is related to the spin 
operator via v = cW/<9k = —vpr x z, which allows us to 
use the TKNN formula for conductivity llj and write 



(15) 



where S 1 ^ is the Kroneckcr delta and we have used the 
fact that the longitudinal conductivity is zero. Hence 
-^cs l = {A/M2d)5eM2v This result is reminiscent of 
the current-induced effective field in single-domain metal- 
lic ferromagnets that belong to the gyrotropic crystal 
class [Hj]. Some significant differences between Ref. [l3| 
and the present work are that Hcs.i (i) docs not depend 
on the strength of spin-orbit interactions in the ferromag- 
net or at the interface (Eq. ([5]) involves vf rather than a 
"spin-orbit velocity"), (ii) reverses sign when il z — > — Q, z 
and vanishes when f2 z = 0, (iii) exerts a dissipationless 
torque provided that the ferromagnet is insulating. 

Hcs, 2 = —{&H/M2D)(A/evF) 2 z x d t & is associated 
with the change in the spin response function under a 



magnetic field (Fig.[2]D): 
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Xm,b{<!) = tLL 



,/*k,n /k+q^ 



k n.n' 



k.n 



(16) 

where q = (lu, q) is the energy /momentum of the magnon 
and r l nn , = (n, k\r l \n', k + q). At q = we get 



(A/ev F ) 2 (-iuj)a H e i: > , where e xy = -e yx = 1 



0. Thus H, 



CS,2 



sim- 



ply increases (if A > 0) or decreases (if A < 0) the Berry 
phase of the isolated ferromagnet thereby renormal- 
izing the parameters entering Eq. |4|. 

When the magnetization of the ferromagnet is uniform 
Eq. fT^]) captures the entire current-induced spin torque 
for weak electric fields. In presence of inhomogeneous 
magnetic textures, one must add the ordinary spin trans- 
fer torque. The microscopic theory for v s - Vf7 amounts to 
evaluating the change of the xy spin-spin response func- 
tion [§] under an electric field (Fig. (2fc). Starting from 
Xm b(q)i perturbing the matrix elements of the spin op- 
erators to first order in E [l5| and expanding the resulting 
expression to first order in q we find (numerically) that 
v s • q oc fl z (E x q x — E y q y ). Furthermore, for realistic pa- 
rameters the torque exerted by Hcs is found to dominate 
over v s • q by an ample margin even when |q| ~ nm" 1 
(note that Hcs does not vary as Ct is slightly tilted away 
from z) . 

Current-induced magnetization switching. — As ex- 
plained above, H CS1 modifies the anisotropy field of 
the ferromagnet in the presence of a Hall current ]h = 
cr H z x E: 



A" 



Q z z ■ 



2D 



ev F M 2 T> 



z xjff, 



(17) 



where K is the anisotropy energy per unit area for the 
magnetic ultrathin film in electric equilibrium. When 
E = the magnetization of the ferromagnet points along 
z. After turning on the electric field, the magnetiza- 
tion begins to preccss around H an and (assisted by the 
damping) equilibrates along the modified easy axis. For 
instance, in a Corbino disk geometry depicted in Fig. [1] 
the electric field produces a crown shaped magnetiza- 
tion. Provided that quantum coherence is preserved, this 
configuration hosts [16| a circulating spin-current propor- 
tional to M(0) xM(</>+<5</>) oc fl z (j H x z) + 0(E 2 ), which 
is radially polarized and persistent (dissipationless). <fi is 
the azimuthal angle around the disk. 

If 3h si evpK/A, O reaches the interface (0 2 = 0) 
in the course of the precession. At that moment, ac- 
cording to Eq. (fT0|) . Ci = and hence d t & = 0; yet 
this is an unstable fixed point and an infinitesimal in- 
plane magnetic field suffices to kick the magnetization 
towards f2 2 < 0. Once this occurs the electric field may 
be turned off and the magnetization will equilibrate to- 
wards — z. Thus a 180° magnetization switching may 
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be completed by combining a dissipationlcss Hall charge- 
current with a very small magnetic field. Nevertheless, 
achieving ju > evpK/A in real materials presents chal- 
lenges. First, jff (E) cannot be larger than ~ lA/m 
(0.5mV/nm) because otherwise the dissipationless quan- 
tum Hall effect will break down [f|. Second, we require 
relatively small coercive fields: H COCT = K/M^d ;$ 0.02T. 
While such a soft perpendicular anisotropy is inadequate 
for the magnetic recording industry, it may find applica- 
tions in magnetic random access memories and magnetic 
field sensors [l8[ . Third, the thickness of the fcrromagnet 
needs to be comparable to the penetration depth of the 
Dirac fermions into the ferromagnetic insulator (< lnm). 
While ultrathin films are commonplace in metallic fcrro- 
magnets [It}, insulating ferromagnets such as EuO or 
EuS present additional experimental difficulties (but see 
Ref. 1 19 for recent progress). Alternatively, one could elec- 
trically manipulate the spin textures caused by mag netic 
impurities placed on the surface of the TI [20, |21| . Us- 
ing A = JM 2 d, 3h = lA/m, Vp = 5 x 10 5 m/s and 
-ffcocr = 0.01T we estimate J > 50meVnm 2 as the condi- 
tion for magnetization switching. Hence J/ a 2 > 0.2cV, 
where a ~ 0.5nm is a typical lattice constant for the topo- 
logical insulator. J /a 2 ~ 0.2eV is an a priori reasonable 
value [13] for the exchange integral between the localized 
moments of the ferromagnetic insulator and the surface 
states of the TI. For stronger perpendicular anisotropies 
(say H cocr 0.05T) the exchange integral would need to 
be of the order of a few eV, and at such strong coupling 
the surface states of the TI would be altered in a way 
not captured by Eq. ([5]). From the precession frequency 
w prcc ~ J-IbH an / h ~ 1GHz we infer switching times of 
the order of a nanosecond. 

There has been some interesting recent work along 
the lines of the above discussion, albeit in topologically 
trivial materials [22j. Two salient differences between 
Ref. (22l | and the present work are: (i) the microscopic 
origin of the change in magnetic anisotropy: in our case it 
is the current-induced spin-polarization of massive Dirac 
fermions (the topological inverse spin-galvanic effect), 
whereas Ref. [22I concentrates on the electrical manipula- 
tion of the atomic positions and distortions of the charge 
distribution; (ii) symmetry of the anisotropy mechanism: 
in our case it is odd under time reversal (because ]h is 
odd), whereas in Ref. [[^ it is even under time reversal 
(because E and charge density are even). 

Conclusions. — When a ferromagnetic film with per- 
pendicular anisotropy is placed on top of a topologi- 
cal insulator, a Hall-current-induced spin torque arises 
which modifies the magnetic easy axis. The origin of 
this torque can be traced to a topological counterpart of 
the long-known inverse spin-galvanic effect, and occurs 
because the electric-field induced magnetization of the 
surface states in the topological insulator is misaligned 
with the equilibrium magnetization of the ferromagnct. 
In Corbino disk geometries this effect might be exploited 



to generate crown-shaped magnetic textures and (for ap- 
propriate material parameters) to switch the magnetiza- 
tion of a ferromagnet by 180° without Joule heating. 
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